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' The Large Electron- Positron (LEP) collider experiments have constrained the mass of the Standard 

Model (SM) Higgs boson to be above 114.4 GeV. This bound applies to all extensions of the SM 
5-H , where the coupling of a Higgs boson to the Z boson and also the Higgs decay profile do not differ 

much from the SM one. However, in scenarios with extended Higgs sectors, this coupling can be 
made very small by a suitable choice of the parameters of the model. In such cases, the lightest 
CP-even Higgs boson mass can in turn be made very small. Such a very light Higgs state, with a 
mass of the order of the Z boson one or even smaller, could have escaped detection at LEP. In this 
work we perform a detailed parton level study on the feasibility of the detection of such a very light 
Higgs particle at the Large Hadron Collider (LHC) in the production process pp — >■ hj — >■ t'^t'J, 
f~| . where j is a resolved jet. We conclude that there are several models where such a Higgs state could 

^ l' be detected at the LHC with early data. 

I) : 

^ . I. INTRODUCTION 

CN 

^ , The Higgs mechanism gives rise to clear and detectable signatures that can be probed at hadron 

' colliders. The Tevatron is now leading a race that the LHC will soon join, to find this spinless 

l/^ , particle. The most relevant searches for a SM Higgs boson by the LEP experiments [11 were based on 

the associated production mechanism, via e~^e~ — >■ Zh{~^ hh) and e^e~ Zh{-^ T~'"r~). A limit of 
\ rrih > 114: A GeV was obtained combining all LEP analyses based on such searches. As the coupling 

psj ■ of a Higgs boson to gauge bosons is fixed in the SM, it is clear that a lighter Higgs can only exist in a 

model where its couplings to gauge bosons (and possibly those to fermions) are reduced relative to the 
On . SM. Moreover, any complementary process, such as e+e^ — > Ah, where A is a pseudo-scalar Higgs boson 

' and h is the lightest Higgs scalar boson (i.e., with the same quantum number as the SM state) in the 

model, which appears, e.g., in a pure 2-Higgs Doublet Model (2HDM) or in the Minimal Supersymmetric 
Standard Model (MSSM), has to be kinematically forbidden. This is a consequence of the obvious sum 
rule g^^zz + 9hAZ valid in a pure 2HDM and in the MSSM. There were also searches at LEP based on 
■ the Yukawa processes e+e~ — 66/i(— > r+r^) in the Higgs mass range ruh =4—12 GeV in Q and in 

the channels bbbb, bbT'^T~ and t~^t~t~^t~ for Higgs masses up to 50 GeV in These studies are only 
relevant for large values of the Yukawa coupling constants and will be discussed in section ITVl 
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Such a scenario, with a very light Higgs, can easily arise by adding Higgs scalar singlets and/or doublets 
to the Higgs sector of the SM. Starting from the first simple extension that allows for a lighter Higgs 
state - i.e., adding a neutral singlet - we will explore all possible scenarios to a maximum of a Democratic 
3-Higgs Doublet Model (3HDM (D)) [5J plus one neutral scalar singlet. We will see that from the very 
simplest extension to all models built thereafter, a very light Higgs state is allowed via a reduction of the 
couplings to gauge bosons. As expected, the more freedom is added by the new fields the more parameter 
space is available to accommodate a light Higgs boson. The possibility of the existence of a very light 
Higgs boson, with a mass below 80 GeV, in the context of models with two Higgs doublets, was discussed 
in [6|. More recently, a very light Higgs boson in the context of the MSSM, with a mass as low as about 
60 GeV, was discussed in 

Assuming this light object would have eluded all LEP searches due to a small enough ghvv coupling 
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(where V = W,Z), the production mechanisms that involve such ghvv couphngs, hke vector boson 
fusion or Higgs-strahlung, have consequently negligible cross sections. Therefore, we have to rely on the 
production modes where Higgs couplings to fermions are present. Under these conditions, the process 
with the largest cross section is gluon fusion. Hence, we have performed a detailed parton level study on 
the feasibility of the detection of a very light Higgs state (below ~ 100 GeV) at the LHC in the production 
process pp — > hj — t~^t~ j Q, where j represents a resolved jet (that indeed proceeds mainly via gluon 
fusion). We have recently presented a similar study for the SM Higgs boson j2] (see also ^0(|) for the 
LHC and a study for the Tevatron was performed in 

The reason to choose such a final state, involving r's, amongst those accessible at hadron colliders, is 
clear. Had we chosen h ^ bb instead, we would have been overwhelmed by QCD background while the 
final state with /i — >■ 77 was not considered because when we vary the Higgs boson mass from 120 GeV 
(where the Branching Ratio (BR) roughly peaks for a SM-like Higgs) down to 20 GeV ^ the corresponding 
BR drops by a factor of ^ 10. Therefore, unless one is exploring a model where the Higgs decay to 
photons is enhanced, the h — > r+r^ mode is the most appropriate channel for light Higgs states. Other 
production channels, like pp — ^ htt, will be explored in the future [l^]- We will show that such a very light 
Higgs could be detected with this process in several extended models and that for particular scenarios 
an early detection is also possible. At the very least, an effort should be made to definitely exclude such 
a light particle and the LHC definitely has the means to do so. 

The plan of the paper is as follows. The next section is devoted to describe the signal and background 
processes in the SM while the following one extrapolates our findings to a variety of beyond the SM 
scenarios with an enlarged Higgs sector. Sect. IV introduces the experimental and theoretical bounds 
enforced in our analysis. Final results are presented in Sect. V (for various scenarios separately) while 
Sect. VI draws our conclusions. Finally, one appendix will help us in classifying four types of 2HDM. 



II. SIGNAL AND BACKGROUNDS IN THE SM 



In a recent work Q we have performed a detailed parton level study on the feasibility of the detection of 
a Higgs boson in the gluon fusion process pp{gg + gq) — hj — t^t~ j at the LHC. In this section we will 
extend the study for Higgs masses below 100 GeV to probe extensions of the SM where a very light Higgs 
is still allowed. The results will be presented considering the case where all Higgs couplings to fermions 
are the SM ones so that they can be used in extensions of the SM with the same final state. As a detailed 
discussion of the parton level analysis was already presented in Q, here we will just highlight the main 
points and refer the reader to reference for details. The signal, pp gg{q) — >■ hg{q) — ^ T'^T~g{q), is a 
one loop process with partonic contributions from gg — >■ hg, gq — ^ hq, which is approximately 20 % of the 
total cross section, and qq — hg, which was shown to be negligible fl3] and was not taken into account 
in our study. We note that this is a parton level study: effects of initial and final state radiation as well 
as hadronisation were not taken into account. 

The SM signal and all background processes were generated with CalcHEP [IJ] and cross checked with 
MadGraph/MadEvent The Higgs BRs to r+r" were evaluated with the HDECAY ^ package (and 
modifications thereof). In the models with an extended scalar sector, the one loop amplitudes for the 
signal pp — ggjq) — » hg{q) — )■ T~^T~g{q) were generated and calculated with the packages FeynArts [l7| 
and FormCalc [18| . The scalar integrals were evaluated with LoopTools [l^ and the CTEQ6L parton 
distribution functions [1^ were used. The jet (leptons) energies were smeared according to the following 
Gaussian distribution 

(1, 

to take into account the respective detector energy resolution effects, where 0.5 is the factor for jets while 
0.15 is the corresponding factor for leptons. 

Each T can either decay leptonically or hadronically. As a three jet final state is very hard to identify 
at a hadron collider, we will concentrate on the other two possibilities - two taus decaying leptonically 
(U) or one tau decaying leptonically and the other hadronically (Ij). These are also the final states with 
robust trigger signatures - the events are selected by an isolated electron with p^ > 22 GeV or an isolated 
muon with > 20 GeV. In both analyses we have considered the main source of irreducible background: 
pp — > Z/j*j — > Uj for U and pp — >■ Z/"f*j Ijj, where one jet originates from a tau, for the Ij case. In 
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pp — )■ Z/"f*j — >■ llj we include all possible combinations of I = e, fj, and in pp — ^ Z/"f*j — >■ Ijj only the 
intermediate state t~^t~ j is included - the jjj signature, where a jet would fake a lepton with a given 
probability, is taken into account in the jjj background. 

The main source of reducible background for the U analysis comes from pp —>■ W~^W~j while for the 
Ij case it is the process pp — >■ Wjj that dominates The tau reconstruction efficiency was taken to 
be 0.3 and accordingly we have used a tau rejection factor against jets as a function of the jet px using 
the values presented in the ATLAS study in [2l|. Finally, we have included the pp tt background 
taken at Next-to-Leading-Order (NLO). By vetoing the events if the tagging jet is consistent with a 6-jet 
hypothesis for \ri\ < 2.5 we were able to discard most of the ti background. The ti background is larger 
for II as there are more possible combinations when the W bosons decay leptonically. 
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FIG. 1: Transverse missing energy distribution for signal and backgrounds for a Higgs mass of 60 GeV. On the 
left for the II analysis and on the right for the Ij case. In this figure all cuts described in the text were applied 
except for the missing energy cut. 



The two other a priori very important sources of reducible background can be brought to a manageable 
level through a judicial cut in the transverse missing energy. This is clearly seen in fig. [1] where the 
transverse missing energy (imbalance of all observed momenta) distribution for a Higgs mass of 60 GeV 
is shown. The huge QCD pp — )> jjj background drops five orders of magnitude as soon as we cross 
the 25 GeV threshold for the missing energy. There is still a tail due to the leptonic decays of c and 
6-quarks which involve a considerable amount of missing energy as is clear from fig. [TJ In accordance 
with CMS [l^l and ATLAS studies ^] we have used 0.001 as the probability of a jet faking one electron 
and, as explained earlier, we have taken the tau reconstruction efficiency to be 0.3 and accordingly we 
have used a tau rejection factor against jets as a function of the jet px [21| that range from 0.01 to 0.001. 

The identification of the Higgs boson signal can only be accomplished by an effective reduction of the 
dominating irreducible Zj background. Therefore, the reconstruction of the mass peak nirT at rrih is 
essential. For a more detailed discussion see H, [ll| . The reconstructed mass distribution is presented 
in fig. [5] for a Higgs mass of 60 GeT^ after all cuts for II on the left panel and for Ij on the right one. In 
both analyses we have sharp mass peaks for the signal and also clear peaks at mz for the Zj background. 

Here we summarise our analysis by specifying the following event selection procedure. 

• We require one electron with p^ > 22 GeV or one muon with pt^ > 20 GeV for triggering purposes. 
An additional lepton in the event has pfp > 15 GeV and pt^ > 10 GeV. A 90 % efficiency is assumed 
for the reconstruction of the electron and muon and the separation between leptons and/or jets was 
chosen as Ai?j(;)j(;) > 0.4 and \rii\ < 3.5 for all leptons. 

• We require that at least one jet has > AOGeV and \T]j \ < 4. 
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FIG. 2: Reconstructed mass mrr distributions for r r decaying leptonically on the left and semi-leptonically 
on the right for a Higgs mass of 60 GeV . 



• We require that the hadronic tau has pij, > 20GeV and |77j| < 4. 

• We veto the event if there is an additional jet with > 20GeV and \rij \ < 5. 

• We apply a mass window rrih — 15 GeV < < + 15 GeV. 

• Events are vetoed if the tagging jet consistent with a 6-jet hypothesis is found with \ri\ < 2.5 (we 
assume a 6-jet tagging efRciency of 60 %). 

• Finally, we require the transverse missing energy to be > 30 GeV. 

In tab. I we present the SM signal and sum of all background cross sections, signal-to-background 
{<^s/^b) ratios and the significance [as/ ^/o^) as a function of the Higgs mass. In the first and second 
columns we show the results for the II analysis while columns three and four are for the Ij case. In 
the last two columns we present the combined values for the signal-to-background ag / cb ratio and the 
significance osl^f^- This study could still be extended to values below 20 GeV, as an experimental 
analysis could well be carried out for such very small values of the Higgs mass. However, the computational 
tools we are using here are not reliable in this Higgs mass regime, so we refrain from investigating this 
phenomenological possibility now. 

This study could still be extended to values below 20 GeV provided the experimental analysis could 
be carried for very small values of the Higgs masses. 

It is clear that the signal observation can be systematically challenging for the larger Higgs masses, 
but the values of the crs/crs ratio can be improved at the expense of the significance by shrinking the 
Higgs mass window especially when its mass is close to the mass of the Z boson. The highest significance 
and the highest as/fJB ratio takes place for mh = 20 GeV because it is the value farthest away from the 
irreducible Zj background. 

In tab. HI] we present the luminosities required for a 95 % Confidence Level (CL) exclusion, 3(7 and 5(7 
discovery of a Higgs boson with SM-like Higgs couplings to the fermions at y/s = 14 TeV as a function 
of the Higgs mass. A light Higgs boson with SM-like couplings to the fermions can be excluded at 95 % 
CL in the mass range 20-60 GeV with less than 1 fb~^ of total integrated luminosity. 



III. EXTENSIONS OF THE HIGGS SECTOR 



Which are the simplest extensions of the Higgs sector of the SM that can accommodate a very light 
Higgs boson? One can write an extensive list of models with very light (pseudo)scalars - it is enough to 
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Mass (GeV) 






aq ffb) 




o-q/o-r (%) 


rro /^//tr (\/fb 1 


20 


11.6 


14.1 


5.4 


28.8 


84.1 


3.24 


30 


11.8 


23.0 


9.9 


35.4 


58.4 


2.97 


40 


11.5 


27.1 


10.6 


41.2 


49.7 


2.76 


50 


11.3 


30.4 


10.9 


47.9 


43.4 


2.58 


60 


11.9 


41.2 


11.3 


63.3 


34.0 


2.34 


70 


13.0 


169.5 


12.0 


149.0 


11.2 


1.41 


80 


13.8 


890.0 


12.9 


856.2 


2.2 


0.64 


90 


UA 


1178.3 


14.1 


1145.6 


1.7 


0.60 


100 


14.9 


1124.7 


15.5 


1142.6 


1.9 


0.64 



TABLE I: Cross sections for signal and sum of all backgrounds after all cuts as a function of the Higgs mass for 
a Higgs boson with SM-like couplings to the fermions. In the first and second columns we show the results for 
the U analysis while columns three and four are for the Ij case. In the last two columns we present the combined 
values for asjoB and as/ ^/(Jb, summed under quadrature. The analysis was done for masses between 20 and 
100 GeV. 



Mass (GeV) 


95 % CL exclusion L (/fo-^) 


3(T discovery L [fb ^) 


5a discovery L [fb ^) 


20 


0.38 


0.86 


2.38 


30 


0.45 


1.02 


2.84 


40 


0.53 


1.18 


3.28 


50 


0.60 


1.35 


3.76 


60 


0.73 


1.64 


4.56 


70 


2.02 


4.56 


12.7 


80 


9.76 


22.0 


61.0 


90 


11.4 


25.7 


71.3 


100 


9.87 


22.2 


62.0 



TABLE II: Integrated luminosities needed to reach a 95 % CL exclusion, 3a and 5a discovery for a Higgs boson 
with SM-like couplings to the fermions, at the LHC. Luminosities shown are for the combined results of the two 
analyses (leptonic and semi-leptonic final states) - signal to background ratio and sensitivities summed under 
quadrature. 



enlarge the parameter space by adding an arbitrary number of fields to accomplish such a goal. However, 
we want these models to reproduce the SM results and to have a high predictive power at the LHC. For 
simplicity we will restrict ourselves to models where Charge and Parity (CP) is conserved in the Higgs 
sector and where natural flavour conservation is assumed 23]. In the SM, the Higgs couplings to gauge 
bosons are fixed by the gauge structure and the Higgs Vacuum Expectation Value (VEV). One of the 
simplest extensions of the SM scalar sector is to add a neutral singlet (i.e., 7 = and Y — 0), where 
/ is the isospin and Y is the hypercharge. Such a singlet does not couple to gauge bosons nor does it 
couple to fermions. The CP-even component of this singlet can mix with the CP-even Higgs field from 
the doublet. Therefore, the coupling to fermions and gauge bosons can only change due to the rotation 
angle related to this mixing. Calling the rotation angle /(%), the couplings are then just redefined as 

#A fix) g'v'A, g'ffi ^ fix) g'ffL (2) 

where V stands for a gauge boson (i.e., V = W, Z) and / is a generic fermion ^. In this case, having a 
light Higgs implies that /(x) ^ 1 and therefore not only the searches based on production and decay 
processes that proceed via couplings with gauge bosons will yield negligible rates but the same is true for 
the ones that rely on fermion-Higgs couplings. Such a light Higgs state h could only be produced then 
in processes involving sclf-Higgs couplings like for example in the gluon fusion or vector boson fusion 



^ If only one singlet is added, /(x) = sinx (or cosx) depending on how one defines the rotation angle. 
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reactions, gg ^ H ^ hh or in qq — )■ qqH — >■ qqhh, respectively. In this scenario, the other CP-even 
Higgs boson, i?, is SM-like in its couplings to vector bosons and to fermions. As for decays, the light 
Higgs BRs are the SM ones because all Higgs couplings fermions and gauge bosons are rescaled by the 
same factor, with the decay to light fermions (6's, c's and r's) dominating by virtue of the small Higgs 
mass that the gauge boson decays are not open yet. Obviously, this scenario cannot be studied with the 
analysis presented in this work because of a negligible production rate for the h state. 

The next step is to add one doublet to the SM to obtain what is known as a 2HDM. Here we will not be 
concerned about any specific 2HDM potential - we just require CP conservation in the Higgs sector. The 
couplings to gauge bosons are universal and we define the coupling of gauge bosons to the lightest Higgs 
as sin(/3 — a) gsM, where (3 is the mixing angle in the CP-odd and charged sectors (tan /3 is also the ratio 
of the Higgs VEVs) and a is the mixing angle in the CP-even Higgs sector. The Yukawa Lagrangian can 
be built in four different and independent ways [l^] if Flavour Changing Neutral Currents (FCNCs) are to 
be avoided. Regarding the Yukawa Lagrangian, there are two clearly different scenarios to be considered. 
The first one is a SM-like scenario where only one doublet, say 02, gives mass to all fermions usually 
referred to as type I model. The second class of models is the one where both doublets participate in the 
mass generation process. With natural flavour conservation, one can build the following models: type H 
is the model where 02 couples to up-type quarks and <j)i couples to down-type quarks and leptons; in a 
type HI model 4>2 couples to up-type quarks and to leptons and 0i couples to down-type quarks; a type 
IV model is instead built such that 02 couples to all quarks and 0i couples to all leptons. We present all 
Yukawa couplings in appendix A. Adding extra neutral singlets to 2HDMs amounts to a redefinition of 
the couplings to the SM particles equivalent to ([2]), that is 

„2HDM , f „2HDM „2HDM , f „2HDM (n\ 

9vVh ~^ Jx^ 9vVh ' 9ffh ^ Jx^ dffh ' \'^) 

where f^- depends now on the number of extra singlets that are added and on the explicit form of the 
scalar potential. 

A class of models which constitute a simple extension of the 2HDM arc the ones obtained by adding 
an arbitrary number, n, (we will call these models 2HDM-|-nD) of doublets that do not couple to the 
fermions. In what follows we will follow closely the discussion in (and notation of) [5^, where a thorough 
analysis of all models discussed in this work is presented. Now we have to distinguish between type I and 
the other types II, HI and IV. In type I models, only one doublet gives mass to the fermions - one can 
then build a new field which is a linear combination of all remaining (n — 1) doublets that do not couple 
to the fermions. If again we choose 02 to give mass to the fermions and 0i as the combined field with 
VEVs V2 and vi , respectively, we will have 

vf + vl = v^uj^, 0<uj<l, (4) 

where w is a function of the remaining (n — 1) VEVs. w = 1 is the 2HDM case while a; < 1 is a signal 
that a non-zero VEV is carried by the linear combination of the (n — 1) fields orthogonal to the light 
Higgs state. 

The situation is slightly more complicated for models II, HI and IV. Taking Model II as an example 
and following [2^] we examine the case where just one more doublet is added, and parametrise the mixing 
with the extra doublet in terms of an angle 6, 

h = cos9h' + sinOho, (5) 

where h' is the usual 2HDM lightest CP-even Higgs boson defined as h' = cos a0,j — sin a4>d in terms 
of the original doublets <&„ and responsible for giving mass to the fermions; /iq is the CP-even state 
of the new doublet 0o that does not participate in the process of mass generation. There is no mixing 
between the two fields when sin 6 — and in this case h — h' . With the usual definition for tan /? = W2/W1 , 
we define cosfi = ■\/ (uj + and sin 17 = vq/v, with < < Tr/2 to write the couplings to gauge 

bosons as 

gvvh = (cos Vl cos 9 sin(/3 — a) + sin Vl sin 6) gyvh (6) 

and to fermions 

cos 6* cos a gj^j COS0 sin a 

gauh — t:—. — o9ffh^ 9ddh — 9iih — p: a9ffh- (') 

cosiz smp cos Si cos p 
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From the physical point of view, adding more doublets will not bring anything new to the particular case 
we are discussing. The same is true if we add an arbitrary number of singlets - the effect is the same as 
for the SM - to reduce all Higgs couplings to the 2HDM fields by the same amount. 

Finally we discuss the case where the fermions masses arise from couplings to three different Higgs 
doublets. We restrict our study to the democratic model described and constrained in [5|, where up- 
type quarks, down- type quarks and charged leptons all get their mass from a different doublet. This 
model is known as the democratic 3HDM and will be represented by 3HDM(D). Following [2^ we define 
tan/3 = Vu/vd and cosO = sinfi = vi/v, where Vd and vi are the VEVs of the doublets 

that couple to the up-quarks, down-quarks and charged leptons, respectively. With these definitions the 
couplings to gauge bosons are 

gvvh — (cos cos 6 sin(/3 — a) + sin Vl sin 6) gyvh (8) 

while the ones to fermions can be written as 

cos 61 cos a gj^.j cosO sin a gj^.j sin 6* ^^.^ 

Quuh — t:—. — adfrhj 9ddh— p: adffhi 9iih — ~ — 7^9ffh- (9) 

cos i I sm p ■' ■' cos S I cos p ' ■' sm \l ■' ^ 

Again, one can add more singlets and doublets to the democratic 3HDM but this will just increase our 
freedom to have a light Higgs boson. We refer the reader to 25] for a discussion on extensions of the 
3HDM. 

In tab. IIHI we present for each model the cross sections for the processes e+e^ — > Zh and pp — > gh 
relative to the respective SM cross sections. In the last column, the BR(ft, r+r^) relative to the SM 
one is shown as a function of the SM BR to up-quarks, down-quarks and charged leptons. It should be 
noted that while the expressions for the cross sections are exact, the ones for the BRs are only truly 
accurate when the Higgs decay to gluons is negligible or when it proceeds mainly through a top loop (like 
in the SM) in which case it can easily be included in the expression using the respective gu coupling. 



Model 




a-(e+e" Zh) 


^{gg gh) 


'BR{h^ T+T-) 


2HDMI 




sin^(/3-a) 


cos"^ a 
sin^ /3 


^ 1 


2HDMI+nD 




uj'^sm^if] - a) 


1 cos'^ a 
uj'-^ sin-^ 


^ 1 


2HDMII 




sin^(/3 - a) 


Flo op 


[1 + (| _ 


2HDMII-HnD 


(cos fl COS ( 


9 sin(/? - a) + sin fl sin Of 


cos^ 17 


[1 + (ff - 1)B7?„]-^ 


2HDMIII 




sin^(/3-a) 


Flo op 


[1 + (| _ 


2HDMIII-|-nD 


(cos Q cos ( 


9 sin(/? - o) + sin f7 sin Of 


cos"^ 
cos^ ft 


\l + C^-l)BRar^ 


2HDMIV 




sin^(/3 - a) 


cos^ a 
sin^ 


[l + {'-^-l){l-BRr)\-' 


2HDMIV+nD 


(cos S7 cos ( 


9 sin(/3 - a) + sin fl sin 


cos^ 6 cos^ a 
cos'^ SI sin-^ 


[l + {f-l){l-BRr)]-' 


3HDM (D) 


(cos Q cos ( 


9 sin(/3 - a) + sin fl sin 


cos^ 6 jji 


[1 + (4 - l)BRu + (4 - l)BRd]-^ 

9i 9, 



TABLE III: Cross sections for e^e~ — >■ Zh and gg — >■ gh relative to the respective SM cross sections for the models 
discussed in the text. In the last column, the BR(/!, — > t^t^) relative to the SM one is shown; in each row, gt 
refers to the model on that particular row. For 2HDM the coupling are presented in the appendix and moreover, 
adding doublets will not alter the gi ratios therefore expressions for 2HDM and 2HDM-|-nD are the same. In the 
case of the 3IIDM, the Yukawa couplings are shown in Eq. ((9|. 



The function Fioop represents the loop contribution which cannot be written as a function of the SM 
cross section. In the SM, the top loop contribution dominates over the bottom loop one by a factor 
{mt/rrii,)^. This is also true in models type I and IV and their extensions. In all other models the factor 
that multiplies the top loop is different from the one that multiplies the bottom loop. We write this 
function symbolically as 

J^^loop — ]—. a'^loop a loop] ■ U'J; 

smp cosp 

In this case we cannot use the SM results and the process has to be recalculated. 
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IV. EXPERIMENTAL AND THEORETICAL BOUNDS 



In this section we present an overview on the bounds of the extensions of the Higgs sector discussed 
in the previous section. We start by noting that we would hke to keep this study as general as possible. 
Hence, we will disregard bounds that require a knowledge about the specific Higgs potential being used. 
This means that only the bounds stemming from couplings to gauge bosons and to fermions will be 
used. We start with the lightest neutral Higgs state. As discussed in the introduction, a SM Higgs 
boson with a mass below 114 GeV was excluded by the LEP experiments. What we are interested in 
this work is to know under what circumstances a light CP-even Higgs boson from a more general model 
could have escaped detection at LEP. All topological searches based on the processes e+e" HiZ and 
e+e" — H1H2, where Hi can be any CP-even Higgs boson and H2 can be either a CP-even or a CP-odd 
Higgs boson were presented in . Considering a scenario where all other neutral bosons are heavy enough 
to have eluded searches based on the various H1H2 combinations in the generic process e+e~ — > H1H2, 
we just need to be concerned with the bounds from e~^e~ — HiZ. Therefore, in this study, the masses of 
the remaining Higgs bosons have no bearing in our analysis. For the pure 2HDMs, the smaller sin(/3 — a) 
is the more the coupling ZZh becomes negligible. Therefore, a very light Higgs could only have escaped 
detection in a parameter region where sin(/3 — a) — )■ 0. Depending on the Yukawa model chosen, there 
are two bounds we have to take into account from direct searches: one that comes from 

a{e+e- ^ Hf^^"^ Z) BRiNf^^"^ ^ bb) _ . 3 BRjHf^^^ ^ 66) 

a{e+e- ^ Hf^Z) BR{Hf^ ^ 66) ~ '"^ " ^' BR{Hf^ ^ 66) ' ^ ' 

and the other one originates from the process 

a(e+e-->g™Z)BR(gi^HDM^r+T-) _ .3 BRjH'f^^^ ^ t+t') 

a{e+e~ ^Hf^Z)BR{Hf^^T+T-) " P> BR{Hf^ ^ t+t~) ' ^ ' 

The reason being that, whatever the values of the 2HDM parameters chosen are, the main decays for 
a Higgs in the mass region below ^ 100 GeV are to 66 and r+r^ if decays to other Higgs are not 
kinematically allowed. The BRs involved depend mainly on the values of a, tan/3 and mh that are 
constrained to obey 

'bR«'((i^M(t+T-)) V 5-"' 'Ve.p 

where j — I, II, HI, IV, BRf/j is the theoretical 2HDM BR and the subscript exp stands for the experimen- 
tally measured value. Note that the limit cannot be shown in the ((/3 — a), rrih) plane because of the angle 
dependence of the Higgs BRs to 66 and to r+r^. It is straightforward to check that when sin(/3 — a) « 0.1 
there is essentially no bound on the lightest Higgs boson mass. To be more precise, in the above mass 
range, a 100% BR to 66 forces sin(/3 - a) > 0.13 while a 100% BR to r+r" implies sin(/3 - a) > 0.16 
regardless of the Higgs mass. Above these values the limit on sin(/3 — a) strongly depends on the value 
of the Higgs mass. Taking as an example Model II, for sin(/3 — a) w 0.2 the limit immediately jumps to 
ruh > 75.6 GeV . We will therefore use sin(/3 — a) = 0.1 as a benchmark value. Note that in the limit 
sin(/3 — a) = we would have a light gaugephobic [1^ Higgs boson which is again not experimentally 
excluded. Because tan /3 is an important parameter in our analysis we note that for a Higgs mass of 40 
GeV for tan/3 = 1 we have 0.57 < sin a < 0.82 while if tan/3 = 30 we have 0.98 < sin a < 1. For larger 
masses the bounds are obviously relaxed but the most important conclusion is that the allowed values of 
sin a are always positive and for large tan /3 we have sin a w 1. It is important to note that these bounds 
do not depend on the total number of (pseudo)scalars in the model under consideration. In the reminder 
of this section we will discuss the limits that do depend on the number of Higgs states. 

The most restrictive bound for a light scalar is the one coming from the muon anomalous magnetic 
moment [g — 2)^ 27]. A verjy detailed account on the subject, including present status of experiment 
and theory, can be found in [2g]. A detailed study of the new physics contributions would force us to 
redo the calculations and subtract the diagrams where a SM Higgs takes part. In what follows we will 
restrict our discussion to the pure 2HDM case - all further extensions of the scalar sector will add more 
freedom to the model, hence they would be much more unconstrained. 

There are two important contributions from extended models to the muon anomalous magnetic mo- 
ment: a one loop contribution, first calculated for the 2HDM in [2^, and the two loop Barr-Zee contribu- 
tion [s^]. The one loop diagram for the light Higgs is proportional to and will therefore have the 
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SM sign which gives a positive contribution to {g — 2)^. This could help cure the present 3cr deviation 
relative to the SM. However, the two loop contribution is proportional to ghfj.+^- 9hbb ^-^d ghfi+fj.- 9hit 
and for SM-like couplings this amounts to a negative addition that, if large, will increase the difference 
between theory and experiment. In fig.[3]we plot Aa^ = {a'^fP — a*^) + a^^^'^^^^ as a function of tan /3 for 
a light Higgs mass of 40 GeV and uia — 200 GeV. We have included both the one loop and the two loop 
contributions and the calculation is presented for the limit a « /3 to comply with the LEP bound on the 
Higgs mass. We have checked that our results for model H agree with the ones presented in [2^ for the 
same limit. For small values of tan /3 the 2HDM contributions can all be safely neglected. As tan /3 grows 
model H makes the discrepancy between theory and experiment to grow. The most interesting scenario 
is the one in model IV - the leptonic 2HDM. In this model, the new contribution moves the theoretical 
calculation closer to the experimental result. This is a very interesting fact, from all 2HDM this is the 
only one that actually helps to cure the problem. Finally, new contributions in models I and HI do not 
vary with tan /?. 

There are other bounds that constrain the pure 2HDMs which deserve a brief comment. Values of 
tan /3 smaller than « 1 are disallowed both by the constraints coming from (the &-jet fraction in 
e+e" Z ^ jets) [3l|, Is^l and from BqBq mixing J53J. Limits from Bs — >■ fJ,~^fJ.~ are not likely to affect 
the models discussed here - it is sufficient to take a large charged Higgs boson mass to avoid these bounds 
in pure 2HDM models . New contributions to the p parameter stemming from Higgs states [s^ have 
to comply with the current limits from precision measurements [36| : \Sp\ < 10^'^. Again for the pure 
2HDMs, there are limiting cases though, related to an underlying custodial symmetry, where the extra 
contributions to 6p vanish. Since we are in the limit of very small sin(/3 — a) there are two limits to 
consider: one is ruh ~ 'niH± and sin(/3 — a) « while the other is m//± — mj^. As we want a light 
CP-even Higgs we choose TO//± ~ friA- Note however that our results only depend on the mass of the 
light Higgs bosons - any change in one of the other Higgs boson masses do not affect our results. 

As mentioned in the introduction, there were also searches at LEP to test the Yukawa couplings [1, |3| 
based on the channels bbhh, bbr^T^ and t^t^t^t^ for Higgs masses up to 50 GcV. In a pure 2HDM 
we are forced to be in a region where a « /?. Hence, the light Higgs coupling to fermions is either 
proportional to tan/3 or to 1/ tan/3. The limits obtained [sl, I3| are for 
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(14) 



which in the most interesting scenarios is reduced to tan ^ BR'^^^^{h — ?> //) - otherwise the data gives 
no useful bounds on the parameters of the 2HDM. The results can be easily applied to model HI and 
model IV because for large tan/3, BR{h bb) « 100% in model HI and BRih r+r") ss 100% in model 
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IV. In any case, even for a Higgs as light as 20 GeV, the obtained bounds are for model III tan/3 < 21 
and for model IV tan/3 < 62. 

The theoretical bounds related to tree level unitarity [s^l and vacuum stability [11] (boundness from 
below) will not influence our results either. Finally we note that although the pure 2IIDMs play a special 
role here because they are protected against c harg e and CP breaking 39] , the same is not true for SHDMs 
or for Higgs models with even more doublets |40| . 



V. RESULTS AND DISCUSSION 
A. 2HDM I to IV 




tan /? 

FIG. 4: Higgs BRs into r+r" relative to the SM one as a function of tan^ in 2HDM Models I to IV with 
sin(^ — a) = 0.1 and = 40 GeV (and also rrih = 80 GeV for model IV). 

With the LEP bounds forcing sin(/3 — a) to be small for a very light Higgs to still be allowed, the BRs 
and cross sections for pure 2HDMs depend almost exclusively on two parameters which we choose to be 
m/i and tan /3. The BR(/i — > r+r^) has a Higgs mass dependence very similar to the SM one. Therefore, 
big differences can only arise in the tan/3 dependence. In fig. 2] we present the Higgs BRs into r+r" 
relative to the SM one as a function of tan/3 in 2HDM Models I to IV with sin(/3 - a) = 0.1. Model I 
has by construction the SM BRs except for exceptional singular points where a given parameter is null^. 
In model II, the down quarks and charged leptons couple to the same doublet and since these are the 
main Higgs decays in the mass region under consideration, the behaviour is again very similar to the SM 
one. For the remaining two models we should distinguish the region of low tan/3 and the region of high 
tan/3. For tan/3 k, 1 and for model HI, there will be an enhancement in BR(/i — >■ r+r") but not in the 
production cross section. When tan/3 ^ 1 is now model IV that sees its BR enhanced and this growth 
depends on the Higgs mass. Therefore the global trend is as follows: in models I and II BR(/i — >■ r+r~) 
is SM-like in all parameter space discussed, in model HI it is enhanced for small tan ji and in model IV it 
grows with tan /3 when compared with the SM value (for other studies on the different Yukawa versions 
of the 2HDM see (H). 



^ The 2HDM fermiophobic Higgs |41.1 originates from model I by setting a = 7r/2. 



11 




FIG. 5: In the left panel we present the ratio between cr(pp — >■ hg)BK{h — )■ r+r^) in model II and the SM 
(t(pp — >■ hsMg) BR(/i5A/ — > T^T~) as a function of ruh and tan/3 — 1 and 30. In the right panel we show the 
corresponding ratio for model IV now for tan/3 = 2 and 3. In both cases we take sin(/3 — q) = 0.1. We also show 
lines of total integrated luminosity 100 p6^^ and 1 fb^^ for model II and 500 p6~^ and 1 fb^^ for model IV. 



The process pp gg{qg) — !• hj proceeds via a quark loop. As the Yukawa couplings are proportional 
to the quark masses, only top and bottom loops give non negligible contributions to the cross section. In 
the SM, the top loop is always the dominant. In our study there are two cases to consider. In models 
type I and IV the light Higgs couples to up and down quarks with the same strength. In this case we can 
write 

a^^^^ipp^gg^h)^'-^ a'^Hpp ^ gg ^ h) (15) 

sm p 

which in turn means that, in the limit sin(/? — a) — > 0, the relation is approximately 

a^^^^ipp^gg^h)^-^ a'^^^pp ^ gg ^ h). (16) 

As we saw earlier, the available constraints force tan/3 to be above 1. Therefore, in models I and IV, 
the 2HDM cross section is the SM cross section for tan /3 = 1 and then drops like tan^ (3 with growing 
tan /?. This means that for tan /3 = 10 the cross section is 100 times smaller and even if the decay 
to r+T" reaches 100 % it will still be 10 times smaller than the corresponding SM cross section. In 
models type II and III the light Higgs couples to the up quarks as cos a/ sin/3 and to the down quarks as 
— sin a/ cos /3. Hence the contribution from each loop depends heavily on the value of tan /3. The larger 
tan /3 is the more the decay to r+r^ becomes negligible in model HI. In model II the width T{h r+r") 
preserves the SM proportionality to T{h — >■ hh) and as the cross section grows for large tan/3, the ratio 
<j{pp — >■ hg) BR(/i — >■ T+T~) to the SM fj{pp — ?> hsivig) BR{hsM — >■ r+r") will also increase. In fig. [5] we 
present this ratio for model II with tan/3 = 1, 30 (left) and for model IV with tan/3 = 2, 3 (right), as a 
function of the Higgs mass and with sin(/3 — a) =0.1. All cross sections in this section were calculated at 
leading order. As stated before, we have chosen the benchmark sin(/3 — a) = 0.1, which is representative 
of a Higgs whose mass is not bounded by the available experimental data. We present luminosity lines 
of Ifb-^ and lOOpb'^ for Model II and 1/6"^ and bOOpb'^ for Model IV. These lines represent the 
integrated luminosity needed to exclude the model at 95 % CL. Obviously, as the Higgs mass approaches 
the Z boson mass, the required luminosity grows. However, there are regions of the parameter space that 
can be probed with less than 100 of integrated luminosity with the LHC working at an energy of 
14 TeV. The regions easily probed are always for the low mass region (below 60 GeV) and large tan (] 
values in Model II and small to moderate values of tan /3 for Model IV. 
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FIG. 6: In both plots the ratio between a{pp — >■ hg) BR(/i — >■ t^t~) and the SM <j{pp — )■ hsAig) BK{hsM -> t'^t~) 
is shown, muhiplied by the factor cos^Q, for rrih ~ 40 GeV and two values of tan /?, 1 (left) and 30 (right) for 
all extensions of the 2HDM in the limit described in the text. We also present the total integrated luminosities 
100 p6~^ and (left) and 100 pb"^ and 1 fb'^ (right). 



B. Beyond 2HDMs 



When we proceed to more general models, the first step is again to make sure the LEP bounds are 
respected. In tab. IIIII we present the Higgs couplings to gauge bosons for models with more than two 
doublets which is written as 

(cos n cos 9 sin(/? - a) + sin sin f (17) 

where there are two new parameters to consider: sin 57 which is a measure of the contribution of the 
new VEV(s) (it can come from the third doublet or from a combination of more than one) and sin 9 
which determines the amount of mixing of the new CP-even field(s) to the lightest CP-even state from 
the 2HDM. Note that, except for the 3HDM case, just the first two doublets give mass to the fermions. 
There are several ways to make this quantity small enough to avoid the LEP bound. We start by 
considering the scenario where almost no mixing occurs between the new doublet and the remaining ones 
and the new VEV is maximal, which can be translated into cosfi ^ 1 and sin 9 <C 1. In this scenario, 
the fermion Yukawa couplings have to be increased to give the fermions the required masses. All cross 
sections, independently of the Yukawa model under study, are now rescaled according to 

„2HDM , 1 „2HDM ^1 q\ 

cr J— a (18) 

cos^ il 

and this means that there is an enhancement already at the level of the cross section. Both a 
and tan^ are now free to vary in all the allowed range, provided theoretical and experimental con- 
straints are fulfilled. In fig. [S] we present the ratio between (j{pp -> hg)BR{h t'^t~) and the SM 
a{pp — >■ hsMg) BR(/isM — >■ T^T^), multiplied by the factor cos^ for rrih — 'iOGeV and two values of 
tan/3, 1 (left) and 30 (right) for all Yukawa extensions of the 2HDM-|-nD. Note that to get the actual 
value of the cross section one needs to multiply it by 1/ cos^ fl and therefore the numbers will always be 
larger than the ones shown in the figures. We also present several values of the total integrated lumi- 
nosity needed to probe the models at 95 % CL. In the left panel we can see that 100 pb^^ are enough to 
constraint a big portion of model III-l-nD while with 2 fb~^ just marginal regions of the models are left 
untested. Between the two vertical lines, the allowed region of the parameter space for the pure 2HDM 
is shown - where a is constrained due to our choice of sin(/3 — a) =0.1 and tan/3 = 1. In the right panel, 
t&nf3 = 30, it is now models II and IV that have the largest cross sections. We show the 95 % CL lines 
for 100pb~^ and 1 fb~^ total integrated luminosity. For large tan/3 only model III-|-nD and regions close 
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FIG. 7: Ratio between a{pp — > hg) BR(/i — >■ t'^t ) and the SM (t(pp — >■ hsAig) BR(/isAf — > T'*'r ) as a function 
of sinS and tan/? = 3 and ruh ~ iOGeV. Models IV+nD and 3HDM(D) are compared. 



to sin q; = in the other models will not be excluded with a few fb ^ of integrated luminosity. The decay 
to tau pairs is negligible in this limit for the 3HDM(D). 

We now consider the scenario sin ^ 1 and a « /? and we will explore, as an example, model IV+nD 
and the democratic 3HDM. Production cross sections are now rescaled as 

which means that they will now be smaller than the corresponding 2HDM cross sections. In fig. [7] we 
show the ratio between a{pp — >■ hg) BR{h — >■ T+r") and the SM (j{pp — ^ hsAig) BR{hsM — >■ t+t~) as a 
function of sin6' for tan/3 = 3 and m/j = 40GeV. For definiteness we take sin 51 = 0.1 and a = /3. The 
prospect of excluding a light Higgs in model IV+nD is good but even with 5fb~^ of integrated luminosity 
only a small portion of the 3HDM(D) will be probed at 95 % CL. Hence, only several years of integrated 
luminosity will allow us to exclude a light Higgs from a 3HDM(D) in this limit. 

There are other scenarios where the LEP bound could be avoided. One is sin^ << 1 and a ~ /?. The 
first relation eliminates the 3HDM(D) light Higgs decay to leptons. For all 2HDM+nD the cross sections 
are again rescaled via the factor 1 / cos^ as compared to the pure 2HDM. Therefore, the results presented 
in fig. [5] will be rescaled with 1/cos^ il to apply for the corresponding 2HDM+nD models. There is also 
the possibility of having sin 51 << 1 and cos6' << 1. Due to the second relation all production cross 
sections will be severely reduced. The only model that has a small chance to be probed in such a limit 
is the 3HDM(D) because the BR to leptons can be of the order 100 %. 

Finally, there is a completely different class of options where none of the previous limits is required. 
By taking 

9vvh — (cos Q cos 6 sin(/3 — a) + sin fl sin 9) gyyi^ — (20) 

and therefore 



sm{i3 - a) ^ ~ta.nnta.ne , (21) 

we do not need to require any special limit to avoid the LEP bound. As an example, if sin(/3 — a) = 
— tan 51 = — tan 9 = — 1 the LEP bound is avoided and the lightest Higgs from 3HDM (D) will have 
SM-like couplings to fermions. Therefore, the SM results shown in table II, can be applied to the 3HDM 
(D). Note that there is no such limit in pure 2HDM cases where a small sin(/3 — a) is required. 
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Type 


I 


II 


III 


IV 




cos a 
sin/3 


sin a. 
cos ^ 


cos a 
sinp 


sin a 

cos /3 


Oidh 


cos a 
sin/3 


sin Q 
cos 


sin a 
cos 


cos ct 
sin 


O-aH 


sin a 
sin^a 


cos a 

cos 


sin a 

sin 


cos a 

cos 


OtdH 


sin a 
sin/3 


cos a 
cos 


cos a 
cos 


sin a 
sin 


Pe 


-cot/3 


tan /5 


-cot/5 


tan /5 


Pd 


-cot/3 


tan^ 


tan/3 


-cot/3 



TABLE IV: Coupling constants for the fermion-scalar interactions. 



VI. CONCLUSIONS 

We have discussed the possibility of finding a very light Higgs boson at the LHC. The existence of such 
a very light CP-cvcn scalar is severely constrained by the LEP bounds which assume the SM coupling for 
the vertex ZZh. If however this coupling is smaller than the SM one, the bound is relaxed and in some 
scenarios the bound on the Higgs mass even ceases to exist. This can be accomplished by the introduction 
of extra Higgs fields. We have seen that the introduction of a neutral singlet is enough to avoid the bounds 
but for this particular process, the production cross section becomes too small. Next, we have introduced 
an arbitrary number of doublets and singlets with the following restrictions: no FCNC are generated at 
tree level and CP is conserved. We have shown that in pure 2HDMs the limit of very small (/3 — a) is 
the only way to avoid the LEP bounds - if one further extends the scalar sector, several combinations 
of different limits in the parameter space lead to the same result. We have also shown that, whatever 
the model is, a few fb~^ of integrated luminosity are enough to probe large portions of the associated 
parameter space. 

With the LHC running and with the search for the Higgs boson on the way, we should ask ourselves 
what to do if we do not find a SM Higgs boson. It seems clear that we should turn our attention to more 
general potentials and in particular to the ones where a light Higgs boson is allowed. However, even if a 
Higgs boson is found and even if looks very much like the SM Higgs boson, we should make sure that we 
did not miss any other (pseudo)scalar particle potentially present in the data. We believe this work is a 
very important contribution to achieve such a goal. 



Appendix A: Yukawa couplings of a 2HDM with flavour conservation 

In this Appendix we present the Feynman rules for the 2HDM Yukawa couplings. Hereafter, the label 
u refers to up-type quarks and neutrinos whilst d to down-type quarks and leptons. Also notice that 
the Goldstone bosons couple just like in the SM, so we do not report their fermionic interactions here. 
Finally, we define 7l = (1 — 75)/2 and "fR = {1 + 75)/2. Using notation already introduced (apart from 
Vij being the Cabibbo-Kobayashi-Maskawa matrix element in the quark sector and equating to 1 in the 
lepton case), one has (see tab. IV) the following Feynman rules: 
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CiGih 2Mw ^^^^ 
UiUih 



djdih 



UiUiH 



UiUiA 
didiA 



_J£_cos_a 
2Mw sin/3 '"""i 



' 2Mw sin p 



'3 
2Mb 



■OLdHrridi 



2M, 



2 Ah 



didiGo 
UidjH^ 



2M;;^^dn^ 

^^|^^eme,(l+75) 

^^l^l^^j [Pdmd, (1 + 75) + cot /3m„. (1 - 75)] 
^/3em..(l-75) 



diUjH TT-M 



7^V* [Pdmd^ (1 - 75) + cot /3m„^. (1 + 75)] 



-^-.(1+75) 



UidjG'^ 



^^V.o [-^d, (1 + 75) + m.. (1 - 75)] 



^y,*[-m..(l-75) + m„,.(l + 75)] 
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